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Abstract. We compare a heuristic count of components of the center 
variety in degree 3 with the equivalent count obtained from known fam- 
ilies. From this comparison we conjecture that more than 100 unknown 
components exist. 



1. Introduction 

In 1885 Poincare asked when the differential equation 

, ^ x+p{x,y) P{x,y) 
y + q{x,y) ' Q{x,y) 

with convergent power series p{x, y) and q{x, y) starting with quadratic 
terms, has stable solutions in the neighborhood of the equilibrium solu- 
tion (x, y) = (0,0). This means that in such a neighborhood the solutions 
of the equivalent plane autonomous system 

x = y + q{x,y) = Q{x,y) 

y = -X- p{x, y) = -P{x, y) 

are closed curves around (0,0). 

Poincare showed that one can iteratively find a formal power series F = 
x"^ + y'^ + h{x, y) + f4{x,y) + . . . such that 




with Sj rational polynomials in the coefficients of P and Q. If all sj vanish, 
and F is convergent then F is a constant of motion, i.e. its gradient field 
satisfies Pdx + Qdy = 0. Since F starts with x'^ + y'^ this shows that close to 
the origin all integral curves are closed and the system is stable. Therefore 
the Sj's are called the focal values of Pdx + Qdy. Often also the notation 
il2j '■= Sj is used, and the rji are called Lyapunov quantities. 

Poincare also showed, that if an analytic constant of motion exists, the focal 
values must vanish. Later Frommer |Fro34j proved that the systems above 
are stable if and only if all focal values vanish even without the assumption 
of convergence of F. (Frommer's proof contains a gap which can ben be 
closed IvWOSQ 
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Unfortunately it is in general impossible to check this condition for a given 
differential equation because there are infinitely many focal values. In the 
case where P and Q are polynomials of degree at most d, the Sj are polyno- 
mials in finitely many unknowns. Hilbert's Basis Theorem then implies that 
the ideal loo = (si, S2, ■ ■ ■) is finitely generated, i.e there exists an integer 
m := m{d) such that 

Sl = S2 = ■ ■ ■ = Sm{d) = =^ Sj = Vj. 

This shows that a finite criterion for stability exists, but due to the indirect 
proof of Hilbert's Basis Theorem no value for m{d) is obtained. In fact even 
today only m(2) = 3 is known. In |vBK09| we prove m(3) > 13 for complex 
centers. 

The proof for m(2) = 3 is conceptually simple: Compute the first 3 focal 
values as polynomials in the coefficients of P and Q under the assumption 
deg(P) = deg{Q) = 2. The 3 polynomials cut out an algebraic variety 
in the space of all differential equations of degree 2. Then decompose, by 
hand or by computer, this variety into its irreducible components. For each 
component prove that all its differential equations have a constant of motion. 

For d = 3 this approach is not feasible because the polynomials Sj are very 
large. They involve 14 variables and are of weighted degree 2j. For example 
the sq can be calculated with our script s6 available at |vBK10b ] and has 
already 95760 terms. The polynomials Sj, j >7 are hard to calculate. Even 
if we would somehow obtain these polynomials, it is extremely difficult to 
decompose the resulting variety into irreducible components. Even I5 = 
(si, . . . , S5) can not be decomposed by c urrent sy stems. So for d = 3 only 



partial results are known, for example |CRZ97 and |Chr05] . In Zol96 



Zolcidek gives a list of 52 families of differential forms known to have a center. 

Our main tool is a statistical method of Schreyer |vBS05j to estimate the 
number of components of the locus Zj where the first i focal values vanish. 
The basic idea is to reduce the equations Sk modulo a prime number p 
and count the number of Fp-rational points of Zi with a tangent spaces of 
fixed codimension. By the Weil Conjectures |Wei49j . which were proved by 
Delinge |Del74| . we know that the fraction of points 

#{¥p rational points on Zi with codimT^^ = c} 



is equal to 



^-^ + higher order terms 



r 



for a disjoint union of r smooth codimension c varieties. If the components 
of Zi are not smooth and disjoint, this number is expected to be smaller. 
More precisely, if the set of singular points has components of codimension 
Cs in the codimension c components of Z, we expect by the same reasoning 
that the fraction of singular points 

/ /^rxN i^i^v rational singular points on codimc components of Zi} 
7p(smg(Zf)) := ^ -^^ 
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is equal to 

• 1 \ C+Cs 



- ) + higher order terms. 

If is smah with respect to p'^" this error does not change the expected 
number '^,p{Zf) significantly. 

Instead of evaluating the S}. at all possible points, we look at a large number 
of random points and obtain an approximate value of ^p{Zf) that can be used 
to estimate r and therefore give an indication of the number of components 
in codimension c. All this is reviewed in Section [3l 

In Section |4] we apply the above method using our implementation of From- 
mers Algorithm. The resulting estimates can be found in Figure [T] 



In Section [5| we analyse Zolfidek's families in detail. We choose random 
points on each family and apply the same statistic as above. Here we find 
that most families are either parametrizing non-reduced components of Z 
or subvarieties of true components. Only 22 families seem to parametrize 
reduced components of Z. Those components can be found in Figure [6j 

Comparing this to our estimate from Section [4] we find that up to codi- 
mension 7 both counts agree. In codim 8 we found heuristic evidence 
for 4 components in Section [4] while in Zolcidek's list we find 5 such com- 
ponents. This apparent contradiction is resolved by showing that two of 
Zoltidek's codimension 8 families {CR4 and CEs) contain the same differen- 
tial forms. For codimension 9, 10 and 11 the heuristic method predicts many 
more reduced components than those that are contained in Zolcidek's lists. 
We therefore conjecture that there are many more components to be discov- 
ered (see Conjecture 5.26). 

The computations for this article were done at the Gauss Laboratory at the 
University of Gottingen. The source code for the Macaulay2 calculations of 
Section [5] is contained in survey2.iii2 using the packages CenterFocus and 
Frommer. These files and the source code for our C-I-+ Implementation of 
Frommers Algorithm can be found at [vBKlOb] . Macaulay2 is available at 

m- 



2. Preliminaries 



If not stated otherwise we work over an algebraically closed field in this 
paper. 

We write the differential equation 

/ ^ P{x,y) 
Q{x,y) 

as P{x, y)dx + Q{x, y)dy = 0. 

Notation 2.1. Furthermore we denote by 

V the 20 dimensional space of degree 3 differential forms Pdx + Qdy. 
W the 14 dimensional subspace of Poincare differential forms 

{x + P2{x, y) + P3{x, y))dx + {y + Q2{x, y) + Qsix, y))dy 
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Definition 2.2. The group 



G ■= Aff2 := 




M V 
1 



detM^^O^ C GL{3) 



with M = ( mil m22 ) ^"^d V = ( ^2 ) is called the affine linear group. G acts 
on the space of differential forms V by affine linear transfomations, i.e. for 



is called the orthogonal group. 0{2) acts on W since it fixes the linear part 
xdx + ydy = ^Dix^ + y^). 

Definition 2.3. A differential form uj = Pdx + Qdy € V has a zero in a if 
P(a) = Q{a) = 0. We say that Pdx + Qdy has a center at o if in addition 
there exist formal power series /x and F centered at a such that //(a) / 
and = fiu. In this case /x is called an integrating factor and F a first 
integral 

Lemma 2.4. //a; /tas a center at a then duj{a) = 0. 

Proof. If cj has a center at a, there exist and F with dF = /xw as above. 
Applying d to this equation we obtain 

= ddF = {dn)oj + niduj) 

Evaluating at a yields 

= {dfi){a)(jj{a) + iJ,{a){duj){a) = iJ,{a){duj){a) 
since ui{a) = 0. Now ij.{a) 7^ by definition, so we obtain (da;) (a) = 0. □ 

Lemma 2.5. // a differential form ui has a center at a then there exists a 
first integral F at a whose Taylor expansion at a 



satisfies Fafl = Fa,i = 0. 

Proof. UJ has a first integral since it has a center at a. Since in the definition 
of first integral only dF appears one can set Fafi = without loss of gen- 
erality. Now Fa^i are homogeneous polynomials of degree i in {x — Qx) and 
(y — ay). Therefore dFi{a) = for all i^l. Now Fx = a{x — ax) + — dy) 
for certain a and j3. We obtain 




The subgroup 




F = Fa,o + Fa,l +Fa,2 + ... 



adx + /3dy = dFi{a) = dF{a) = {lMo){a) = 



and conclude Fi = 0. 



□ 
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^^'(pf^"''l rK^"' nl=E«.(^'^)(- 
.3 = 1 



Definition 2.6. In the situation of Lemma 2.5, -F2,a =: ^2('^)0) is called 
the quadric associatied to oj in v. The rank of F2{u},a) is invariant under 
affine coordinate tranformations. 

Remark 2.7. If uj has a center at (0,0) and rankF2(w, (0, 0)) = 2 we 
can assume that F has no constant or linear terms as above. Over an 
algebraically closed field we can find a coordinate change such that 

and UJ = xdx + ydy + . . . , i.e u; is a Poincare differential form. 

Over an arbitrary field this is only possible if additional conditions are satis- 
fied. For example over M one must assume that the quadratic form associated 
to F2 is positive definite. 

Definition 2.8. Let Pdx + Qdy be a Poincare differential form of degree 
3 over a field of characteristic 0. One can then use Frommer's algorithm to 
find a formal power series F £ K[[x,y]] with 

00 

2j+2 ^y2j+2y 

In this situation Sj{P, Q) is called the jth focal value of Pdx + Qdy. From- 
mer's algorithm also implies that sj is polynomial on W and has rational 
coefficients. We call Sj G Q\pij, Qij] the jth focal polynomial. 

Remark 2.9. By analysing Frommer's Algorithm |vBC07j one can show no 
prime factor of that the denominator of sj is bigger than 2j + 2. Therefore 
Sj mod p is well defined for j < {p — 3)/2. 

Definition 2.10. We define the ideals 

Ij = {Sl,. . . ,Sj), loo = (Sl, S2, . . . ) 

and their vanishing sets Zj = V{Ij) C W . is a variety whose points are 
exactly the Poincare differential forms with a center at (0,0). We therefore 
call it the center variety. 

Remark 2.11. In the case of degree 3 differential forms considered here, 
Q\pij,qij] has 14 variables. Hilbert's Nullstellensatz implies that loo can be 
generated by finitely many elements, therefore there exist a number m := 
m(3) such that Z^o = Z^ and Z^o 7^ Zm-i- The precise value of m(3) 
is unknown. In [vBK09| the inequality m(3) > 13 is proven for complex 
centers. Since one can not study Z^o explicitly we analyze Zi^ in this paper. 
If m(3) = 13 this is equivalent to analyzing Z^q. Otherwise we have Z^o C 

■2^13- 



3. Finite Field Heuristics 



In this section we explain how one can obtain heuristic information about 
a variety X C A" by evaluating its defining equations at random points. 
For an extended discussion about this method see |vBS05j or jvBOSij . An 
application of this method to the Poincare center problems in some solved 
and some unsolved cases is described in |vB07] . 
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Definition 3.1. Let X C A"'(Fp) be an algebraic variety. Denote the num- 
ber of Fp-rational points of X by |X(Fp)|. Then 

|An(F,)| 

is called the fraction of¥p-rational points of X in A". 

Remark 3.2. If X has r irreducible reduced smooth components of codi- 

mension c and all other irreducible components have larger codimension 

then the Weil-Conjectures imply that 

/l\c 1 
7p(X) = + higher order terms in - 

We will estimate "jpiX) statistically by evaluating the equations defining X 
in a number of randomly chosen points. 

Definition 3.3. Let X C A"(Fp) be an algebraic variety. For a sequence 
S = {xi, . . . ,xn) of Fp-rational points in A"(Fp) we call 

nfp{x,s) — 

the empirical fraction of ¥p-rational points. 

Remark 3.4. The distribution of %{X, S) on the set of all sequences S of 
length N is binomial with mean fi{^p{X, S)) = jpiX) and standard devia- 
tion 



a{%{X,S)) = 



jpiX){l-^piX)) ^IjpiX) 



N \l N 



This allows us to obtain an estimate of 7p(X) and then of r and c by eval- 
uating the equations of X in many random points. More information is 
obtained, if we also calculate the tangent space of X in these random points: 



Definition 3.5. Let X C A" be an algebraic variety defined by /i = • • • = 
fr = 0. Then the tangent space of X in a point x & X is defined as 

Cj ' ' / i=l...rj=l...n 

Remark 3.6. Let X' C X G A" be an irreducible component, x E X' a 
point and Tx',x the tangent space of X' in x. Then 

codimX' > codimrx',a; 

with equality for general points if X' is reduced. We therefore consider only 
points with codim Tx'^x = c in estimating the number of components of codi- 
mension c. By the inequality above we disregard all points on components 
of codimension greater then c. 

These arguments lead us to 

Heuristic 3.7. Evaluate the equations of X in A'^ random points Xi over Fp 
and calculate the tangent spaces Tx,Xi in these points. Then estimate 

^{reduced codim c components} ^ #{^ I codimTx,x^ c}^c 
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with an estimated error 

codimrx,a-^ = c} ^ 

In this paper we have used <I> = 2 to obtain a confidence level of approxi- 
mately 95%. 

Caution 3.8. Let be the subvariety of X whose points have a tangent 
space of codimension c. Then above heuristic means that statistically the 
hypothesis 7p(X'^) = r{l/py can not be rejected with confidence of more 
than 4.6%. Algebraically this proves nothing, but gives a way to arrive at a 
reasonable conjecture about X. 

Caution 3.9. It is possible that X contains a component Y that is 
irreducible over Q but decomposes into several irreducible components 
Yi,...,Yk over the algebraic closure Q, i.e. the Galois group Gal(Q/Q) 
acts transitively on the Yi. Over a finite field a Yi is rational if the Frobe- 
nius endomorphism fixes Yi. The expected number of such components is 1. 
Therefore our heurisic is an indication of the number of reduced irreducible 
components over Q and not over Q or C. 

Remark 3.10. If the components of X are not smooth and disjoint, then 

#{i I codim Tx,x, = c} ^ 
N ^ 

is expected to be smaller than the actual number of reduced codim c com- 
ponent. More precisely, if the set of singular points has components of 
codimension Cg in the codimension c components of X, we expect by the 
same reasoning that the number of singular points on codim c components 
to be approximately 

r,N 



If is small compared to p'^" our Heuristic 3.7 is therefore also useful in the 
presence of singularities. If not, the number calculated can still be used as 
a heuristic lower bound on the number of reduced components. 



4. Experiments 



Using the heuristics described in Section [3] one can estimate the number 
and codimension of reduced components of the center variety Zqo. For this 
we study Z13 ^ -^00 as an approximation. This is possible because From- 
mer's algorithm |Fro34j . |vB07j . |MorOOj provides a fast way to calculate the 
focal values of a given Poincare differential form even though the explicit 
polynomial expressions for the focal values are not known. 

Experiment 4.1. We examined 402376372880300032 = 4.02 x 10^^ points 
over F29 and determined the rank of the Jacobi matrix if the first 13 focal 
values vanished using our implementation of Frommers algorithm |vBKlQb] . 

This would take about 11 years of CPU time on a 2.3 GHz AMD Opteron 
Prozessor with 128 KB Ll-Cache and 512 KB L2-Cache using our newes im- 
plementation of Frommers algorithm and a parametrization for the solution 
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set of the first three focal values to speed up the process. We distributed 
the work to 56 processors. 

The heuristic estimate derived from this experiment is shown in Figure [TJ 
Interesting differential forms found in this and other computer experiments 
as well as statistics about theses experiments are collected in our online 
database j vBKlOa] , 

Remark 4.2. To test our implementation we have used it to recalculate 
the focal values of the examples in [HoOlj . Also the focal values of our ex- 
ample in [vBK09] were calculated independently by Colin Christopher using 
Reduce and agree with ours modulo 29. Furthermore the fact that for most 
Zoladek differential forms we indeed find points whose first 13 focal values 
vanish (see Section [5]) can be interpreted as another test of our implemen- 
tation. 

To test the parametrization of the first thee focal values we compare the 
results obtained with and without using parametrization. 

To ensure that our experiments can be repeated we use a pseudo random 
number generator and store the svn revision number of the program version 
used to do the calculation in our database. 

Remark 4.3. By applying elements of the group 0(2) to a given differential 
form oj over F29 we obtain further differential forms that have exactly the 
same properties as oj. Now the group 0(2) has 2 • 28^ elements over F29 and 
therefore only approximately 

29I4 17 

^ ^ 1.9 X 10^^ 

2 • 282 

fundamentally different differential forms exist over F29. Since we choose our 
points randomly it can happen, that some forms that are equivalent with 
respect to 0(2) have been analysed several times. This makes no difference 
for our statistics, but prevents us from looking at all points even though we 
have made more than 1.9 x 10^'' calculations. More precisely the propability 
of missing a general 0(2) orbit was 

1 X 4.0x10" , 4 (|^i„17x 

for our experiment. Therefore one can expect that we have seen about 
of the fundamentally different differential forms. 



5. Zoladek's Lists 



In [Zol94| and |Zol96| Zoladek has given a list of 52 families 

(/.: A" ^ y 

of degree 3 differential forms with a center. They are divided into 17 rational 
reversible systems and 35 Darboux integrable systems but this distinction 
is not needed for our survey. No claim on the completeness of this list is 
made. 
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Figure 1. Number of Poincare differential forms whose first 
13 focal values vanish over F29 found after shifting through 
402376372880300032 random differential forms. 

Remark 5.1. In |vB07j we proved that Zoltidek's families CR5 and CRi are 
subfamilies of C-D4 and similarily CR12 and CRiq are subfamilies of CD2- 
We will therefore not consider them in this paper. 

Remark 5.2. Notice that the trivial Hamiltonian component of differential 
forms oj that satisfy 

uj = dF 

for a polynomial F of degree 4, is not on Zoladek'slist. 

Remark 5.3. In printing long lists of polynomials it is impossible not to 
introduce missprints. For this paper we have started from the implementa- 
tion in Ulrich Rheins thesis [RheOSj and made further corrections. Together 
we have made the following changes 

• In CRi we changed the second occurrence of (7 to a new variable, 
enlarging the family to all symmetric forms. 

• For Ci?3 we changed (2a^ — 6) to (2a — 6^) in the expression for F and 
5a/2 to 3a in the expression for H. The first change is in |Rhe08] 
the second isn't. 

• For C-R5, X and y have to be exchanged in [ZolQSj. Also / has to be 
changed to ly. This is already corrected in [RheOSj . 

• In CRii a sign mistake was introduced in [RhcOS] 

• For CRn the derivatives r]x and r]y were not calculated correctly in 
tRhe08]. 

• For CDiY the equation 



A(3{(3 - l)a^ + 4/3(3 - 2/3)a + (3 - 2/3)(l - 2/3) = 
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must be satisfied. Fortunately the curve defined by this equation is 
rational and can be parametrized by 

" ~ -2tlti + st^tl - utotl + mf 

-2tlti + Stltf - Utot'l + I2t\ ■ 

We substituted this parameterization into the expression for CDn 
set ti = 1 and considered only the numerator of the resulting ex- 
pression. The parameterization was kindly computed for us by Janko 
B5hm [BohlOj . 

• For C-D24 we did not find any centers over F29 

• In CD25 the coefficient of x'^ was changed from a to a. This misprint 
was already corrected in |Rhe08] 

• In CD26 the division /2 must be erased. This was also found by 
[RheOS] , 

• For CD32 we did not find any centers over F29 

• From CD^3 we obtain degree 4 differentials for generic coefficients. 
Only in the case a = 1 we were able to factor out another factor x. 
We therefore only use CD^z with this additional restriction. 

• Some families can be trivially enlarged by scaling with a nonzero 
scalar. We did this for all CD^s except CD^ and C-Dg by multiplying 
the formula given by Zolcidekwith the variable aaig- 

The families we used are contained in our Macaulay2 package CenterFocus 
IvBKinbj . where we have renamed the variables a, . . . , t, a, 7 to 
aai, . . . , aaig. 

To estimate what part of our statistic in Figure [l] is explained by 
Zolcidek's examples we need to take into account, that Zolcidek's examples are 
general degree 3 differential forms in V while we are interested in Poincare 
differential forms in W. Over an algebraically closed field every degree 3 
differential form uj with a non degenerate center can transformed into a 
Poincare differential form by an affine transformation. It is the purpose of 
this section to formalize this process and keep track of the dimensions of the 
families involved. 

Definition 5.4. The affine linear group G acts on the center variety. There- 
fore if 

(/.: A" ^ y 

is a family of differential forms with a center, then 

tp: GxA" 

ig,a) ^ {g{(t>{a))) 

is a (possibly larger) family of differential forms with a center that is invari- 
ant under action of G. Furthermore 

ImV' n 

is a variety of Poincare differential forms. 
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Figure 2. Dimension of Imcp for Zolcidek's families 



Remark 5.5. Im^ n W can have several components Wi of which at least 
one contains differential forms with a center at (0,0). The subset of such 
differential forms is then dense inside this component. 

Lemma 5.6. Let <j): A"' V be a morphism, and D(p its differential. If 
a E A is a integral point with rank(D(/))(a))Fp = n then dimcj) = n. 

Proof. We have the following inequalities 

n > iai\k{D(j))[a) > rank((D(/>)(a))iFp = n. 

Since the D(f) drops rank only on Zariski closed subsets of A"' we also know 
that D(j) has generically rank n. If follows that (j) is generically locally injec- 
tive and therefore dimlm^ = n. □ 

Calculation 5.7. We compared the number of variables n involved in the 
definition of Zol§,dek's families with the rank of D<j) in a random point a 
using our script rankDif f erential. For all families both numbers agreed. 
Figure [2] contains Zol§dek's families sorted by n = rank Dcp = dim Im (j). 

For the remaining calculations we need the following theorem on the dimen- 
sion of fibers of a morphim: 

Theorem 5.8. Let (p: X ^ Y he a morphism of irreducible varieties over 
an algebraically closed field. Then 

dim. (f)~^{y) > dimX — dimy 

for all y gY and there exist a Zariski open subset U CY such that 

dim(/)~^(u) = dimX — diml" 

for all u G U. In this situation dimX — dimy is called the generic fiber 
dimension. 

Proof jMum88L 8, Theorems 2+3] □ 



Definition 5.9. For a family (p we denote by = n — dim Im (p the generic 
fiber dimension of (j). For all Zoladek families considered in this paper we 
have seen di = in Calculation 5.7 
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Definition 5.10. Let (p: A" V he a family of differential forms and 
a E A" a point. Then 



Ga = {geG\g{(P{a))elm^} 

is called the set of irrelevant elements of G with respect to a. Consider now 
the variety 

X = {(5,a)|5(0(a)) G J^} C G x A" 

and the projection 

vr: X ^ A" 

then 7r^^(a) = Ga- From Theorem 5.8 we obtain that for almost all a 



dimGa = dimX — n =: d2- 
we call ^2 the generic dimension of Ga- 

Calculation 5.11. In Figure |3] we list subsets Ha C Ga for almost all a for 
all of Zoladek rationally reversible families. That these are indeed subsets is 
checked by our script islrrelevant. 

Calculation 5.12. For every rational reversible Zolcidek family 
we calculate dimGao for a random element ao using our script 
ideallrrelevantElementsRandom. The results can also be found in 
Figure |3j Since dim Ga,, is always bigger then the generic dimension of Ga 
we obtain for almost all a G A": 

dim Gao > (^2 > dim Ha 

From the Figure [s] we see that for every Zoladek family these inequalities 
have to be equalities and we can calclate d2- 

Calculation 5.13. For every Darboux integrable Zoladek family 
we calculate dimGa^ for a random element using our script 
ideallrrelevantElementsRandom. We obtain that this dimension is 
zero for all CDi. Since 

dim Gao ^ > 
We obtain (i2 = for these cases. 

Proposition 5.14. Consider 'ip: G x A" -^V as above. Then 

dim Im ^/' = n + 6 — — d2- 

Proof. Let ^{g,a) = g{(j){a)) be a generic element of liaip. Then the fibere 
over this element is 

F = ^-^g{m) 
= {{hMKm) = 9{<i^{a)} 

= {{hM{m = h'^9{m)] 
= {(/i,6)|0(6) = /i(,/.(a))}cGaX A". 
Now consider the projection 



TT-.F^Ga 



A SURVEY OF THE POINCARE CENTER PROBLEM IN DEGREE 3 



13 



Familiy 


H 


dim H 


dim Ga 
(5.12) 


n 
(5.7) 


codim 
ImV' n 
(5.19) 


CRi 


( mil \ 
V m22) 


U) 


3 


3 


10 


6 






^m22 \ 
^ m22y 


(0) 




1 


1 


6 


8 


Gits 


( 


^mii \ 
^ imii) 




1 


1 


5 


9 


O JX4 


f mil 

V 


V2miiaa^ 
V2miiaa^^ + 


mil) \V2j 


2 


2 


7 


8 


CR5 


subfamily of CD^ |vBn7j 


CRq 


/ mil V2miiaa^ 
\ m22 


)( 




3 


3 


8 


8 


CR7 


subfamily of CDi [vB07] 


CRs 




^m22 \ 
V m22y 


(0) 




1 


1 


5 


9 


CRg 




^77122 \ 
^ 77122/ 


(0) 




1 


1 


5 


9 


CRio 




^77122 \ 
V 77722/ 


(0) 




1 


1 


5 


9 




f mil 
V 


W277iiiaa3 
V2miiaa^^ + 


mil) \V2) 


2 


2 


8 


7 


CR12 


subfamily of CD2 |vBn7j 


CR13 




^77722 \ 
V 77722/ 


(0) 




1 


1 


5 


9 


CRi4 




^77722 \ 
^ 77722/ 


(0) 




1 


1 


6 


8 


CR15 




^77722 \ 
V 77722/ 


(0) 




1 


1 


5 


9 


CRw 


subfamily of CD2 [vB07J 


CRn 


(i ?)(o) 








2 


11 



Figure 3. Calculating the codimension of Imip n W for 
Zolcidek's rationally reverible systems 



For heG a we obtain 

TT-Hh) = {b\(t>{b) = hi4>{a))} 
= {b\cf>{b) = </.(a')} 

For generic a and h we therefore have 

dim F = dim Ga + dim (/)~^ (a) = di + d2 
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□ 



by Theorem 5.8, Using Theorem |5.8| again for ip and generic F we get 
dimlmV' = dim(G x A") — dimF = n + 6 — di — d2 

Proposition 5.15. Consider the variety 

X = {{g,uj)\g{co) € W} C G x Im^ 
and Xq G X an irreducible component. Let 

vr : Xq — )• G 

(g,^) ^ g 

be the natural projection. In this situation all fibers of vr are isomorphic 
and iT~^{id) is an irreducible component of Imip n W. Furthermore this 
component has the dimension dimXo — 6. 

Proof. G operates on X via h{g,uj) = {gh, h^^oj). Since G is irreducible 
it also acts on every component Xq C X. With this operation 7r^^(/i) = 
h^^{TT^^{[d)). This proves the first claim. Now 

^-i(id) = {lo\uo G W}nXo. 



This proves the second claim. The third claim follows from Theorem 5.8 □ 



Lemma 5.16. Let X be the variety considerend in Proposition \5.15\ and the 
natural morphism 

rj : X ^ X limp 

Then 

Im?7 = {{v, uj)\uj(v) = duj{v) = and iankF2{v, uj) = 2}. 

Proof If oj' := i^(l\){u}) lies in W, it satisfies uj'{0) = du}'{0) = and 
rankF2(0,a;') = 2. But then uj{-M-^v) = duj{-M~^v) = 0. With v := 
—M^^v this shows 

Imrj C {{v,uj)\uj{v) = du}{v) = Oand rank F2{v , oj) = 2}. 

Conversely consider {v,lo) with uj{v) = duj{v) = and rankF2(?),a;) = 2. 
Then with g = (q we have uj' = g{ui) satisfiying w'(0) = 0. This shows 
that oj' is of the form 

/ 1,1, 

UJ = UJi + (jJ2 + . . . 

with Lo'^ = dF2. Since ranki<2 = 2 there exists an element h = \ ) such 
that /i(F^) is \{x^ + y"^). It follows that 

bj" := h{uj') = xdx + ydy + . . . 

and {h o g,uj) is an element of X with image {v,uj). □ 
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Figure 4. Estimated number of center variety components 
parametrized by Zoiadek's rational reversible families 



Lemma 5.17. In the situation of Lemma 5.16 the fibers of 

rj: X ^Imr] <Z x Im 

are isomorphic (as varieties) to 0(2) C G. In particular the components 
Xi of X are in 1 : 1 correspondence with the components Ei of Im r] and 
dimXj = dimi?j + 1. 

Proof. The group G acts on x Im^ via 

h{v,uj) = {h{v),h-\uj)) 
where for ^ = ( ^' 'l' ) we set 

h{v) := -{M')-\v' + v). 

With this action the morphism rj is G covariant. It follows that the fiber 
r]~^{v,uj) is isomorphic to a fiber rj~^{0,uj') with uj' = h^^{uj) for an h with 
h{v) = 0. We have 

rj-'^{0,uj') = {M\{^^ \ ){uj') e W A -M~^v = 0} 

= {M\{^^^,){c.')eW}. 

This set is non empty, since (0, oj') is in the image of t]. Therefore there 
exists an h' such that h'{0,u') = (0,cj") with u" G W. Now 

V-H0,u^") = {M\r^1)iu")GW} 
= 0(2) c G 

since only elements with MM'^ = ( q 5 ) fix the linear part = xdx + 
ydy. □ 

Corollary 5.18. If (f) is a family of differential forms with a center whose 
generic element has only finitely many zeros and a center of rank 2, then 

dim W n Im tp = n + 1 — di — d2. 
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Figure 5. Estimated number of center variety components 
parametrized by Zoladek's Darboux integrable families 

Proof. By assumption a Zariski open subset of lm(j) contains differential 
forms with a rank 2 center. Since this fact is invariant under action of G 
the same is true for Imip. For every such element u G Imip one can find an 
element g £ G such that g{u}) is a Poincare differential form in W. Therefore 
^ o ?7 is dominant. If lo has only finitely many zeros then ^"^(o;) is finite, so 
^ is generically finite by our assumptions. We obtain 

dim X = dim Im ip + 1. 

Using Proposition 5.14 and Proposition 5.15| we obtain 

dim Im ip DW = dim X — 6 = dim Im ip — 5 = n + 1 — di — d2. 
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□ 



Calculation 5.19. Using Corollary |5.18 we calculate dimVF n Imi/^ for 
Zolcidek's families of rationally revesible centers. The results can also be 
found in Figure [sj For Zolcidek's families of Darboux centers we have di = 
d2 = and therefore the dimensions are equal to n + 1 and can be read from 
Figure [2} 

Remark 5.20. We collect the previous definitions, lemmata and propsitions 
in the following diagram: 



G 



codim 6 Fibers 



Xi 



c X c G X lm^p 



0(2)-Fibcrs 



Im?7 



finite 



Im'i/' 



Im^ 



where labels of the ^- and rj arrows denote the expected fibers. Special fibers 
could have a different structure. 

To estimate the component structure of 7r^-'^(id) = Im^ n we use again 
our heuristic approach. 

Calculation 5.21. Starting from a family (p: A" — )• V, we find rational 
points on Im-i/; DW as follows. First choose a rational point a G A" and 
consider the differential form u = ip{id, a) = {(pia)) E l^^l^p. li vi, . . . ,Vk are 
the rational symmetric zeros of u, then the rational points in the preimage 
of ^ are 

= {{vi,uj) \ i = l,...,k}. 
For each pair {vi,uj) € A^ x Im^ the preimage of t] is 

rj-\vi,uj) = {iM,-Mvi,ujg\iM, -M'^Vi)iio) G W} 

Now 

y 1 V ~ V 1 1 ^ 
and := ( Q ){^) has zeros at {vj—Vi) in particular one at zero. Therefore 
Ui is of the form 

oji = liixdx + li2{xdy + ydx) + l22ydy + higher order terms 

We have {'^^ ^){u;i) £ W if and only if 

M%M = (iO) 

with 

* ^ hi '22 ' 



Since by Lemma 5.17 the solution is a one dimensional space, we can fix 



one entry of M and and generically obtain finitely many rational solutions 
Mil, ■ ■ ■ , Mil. Ill this manner we have found finitely many rational points 

{g,uj) e (r/o^)-i(a;) C X C G x Imcf) 
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Figure 6. Zoltidek families that heuristically parametrize 
reduced components of the center variety 

with 

g £ {{Mij,-MijVi)}. 

To obtain points in H Imip = 7r~^(id) we operate with g^^ on the whole 
situation, and obtain points 

(id,g(a;)) G X C G x ImV 

The points g{uj) can then be analysed with our implementation of From- 
mer's algorithm. If the first 13 focal values of g{uj) vanish we calculate the 
codimension of the tangent space to Zi^ D Z^o in these points. The results 
of doing this for 2000 random choices of a in each of Zol§,dek's families are 
available as hash tables experimentsCR and experiment sCD in survey2.m2. 
For the families CD24 and CD32 we did not find any differential forms this 
way. 

We now want to identify those Zoltidek -families that define reduced compo- 
nents of the center variety Z13 D Z^o ■ For this we use again our finite field 
heuristic. 

Calculation 5.22. Consider a family 0: A" — )• V and set d = Wnlmtp. W 
might have several components Wi of which at least one has dimension d. By 
the procedure above we expect to find approxemately 2000/p'^~'^'™^* points 
on a component Wi of WCilmip. The generic codimension of a tangent space 
to the center variety in points of Wi is 14 — dim Wi if and only if Wi is also 
a reduced component of the center variety. We can therefore heuristically 
identify reduced components of the center variety Z13 by scaling our point 
counts by p'^~^^"'"'^/2000 where c is the codimension of the tangent space a 
each point. The result is contained in Figures |4] and [5] 

Remark 5.23. A family (f> will have all numbers calculated above close to 
zero if one of the following holds 

(1) (j) defines only a subfamily of a true component of the center variety 
and the codimension of the family inside the component is at least 
one. 

(2) (j) defines a non reduced component of the center variety 

(3) The generic point of (p does not have a symmetric center. 
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• p=31, 14 equations, 2.3 x 10'^17 points 
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Figure 7. Our heuristic predicts that up to codimension 
8 all reduced components of the center variety are known. 
For higher codimension many components are waiting to be 
discovered. 



We suspect that all three possibilities actually occur. The third case can 
be easily detected by analysing a generic point. This shows that families 
CD^s and are of this kind. Probab ly t his is e ither due to misprints 

introduced by us or by misprints in (Zol94 or Zol96 that we were not able 
to find and correct. 



To distinguish between the cases (1) and (2) is much more difficult. 



Remark 5.24. Notice that only smooth points on each components have 
the correct tangent dimenesion. Therefore we expect the results of the above 
scaling to be less than 1 for each component of the center variety. We have 
collected those families that do parametrize a reduced component of the 
center variety by this heuristic in Figure [6} 

Comparing the number of components contained in Figure [6] with those of 
Figure [T] we find that up to codim 7 both counts agree in codim 8 there 
are 5 components given by Zoladek , while we see only 4 in our heuristic. 
Fortunately Ulrich Rhein has found numerical evidence for Ci?4 C CRq in 
his Diploma Thesis |Rhe08| . It is not difficult to prove that this is indeed 
the case: 

Proposition 5.25. All differentials parametrized by Zolqdek's family CR4 
are also contained in Zoig,dek's family CRq. 
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Proof. One can obtain Ci?4 from CRq by setting k = and renaming the 
variables as follows r— T-f;— >-p— >-n— p-Zcin Zolcidek's notation. □ 

With this correction we have compared our heuristic component count with 
the components detected among Zol§,dek's list in Figure [t) We observe, that 
up to codimension 8 both counts agree. Starting from codimension 9 there 
seem to exist many more reduced components than previously known. We 
therefore 

Conjecture 5.26. The number of reduced components of the center variety 
in degree 3 is 

• 1 in codimension 5 

• 2 in codimension 6 

• 4 in codimension 7 

• m codimension 8 

• at least 12 in codimension 9 

• at least 33 in codimension 10 

• at least IJ^ in codimension 11 

• possibly further components in codimension 12 

Remark 5.27. The Macaulay2 calculations made in this section are con- 
tained in the file survey2 .m2 using the packages CenterFocus and Fronmer. 
All three are available at |vBK10b] . Macaulay2 is available at |GS| . 
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